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Abstract We present a determination of the distance modulus of the globular cluster NGC 5904 (M 5), obtained by means of 
the subdwarf main-sequence fitting on the (V,V-I) color-magnitude diagram. The subdwarf sample has been selected from the 
HIPPARCOS catalog in a metallicity range homogeneous with the cluster ([Fe/H] m -1.1). Both the cluster and the subdwarfs 
have been observed with the same telescope+instrument+filters setup (namely, ESO-NTT equipped with the SUSI2 camera), 
in order to preserve homogeneity and reduce systematic uncertainties. A set of archival HST data has then been used to obtain a 
deep and precise ridge line. These have been accurately calibrated in the ground photometric system by using the NTT data and 
used to fit the cluster distance modulus. By adopting the most commonly accepted values for the reddening, E(B-V) = 0.035 
and 0.03, we obtain respectively /j = 14.44 ± 0.09 ± 0.07 and /j = 14.4 1 + 0.09 + 0.07, in agreement with recent determinations. 
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1. Introduction 

The ages of the Galactic Globular Clusters (hereafter GGCs) 
represent one of the milestones on the road of our comprehen- 
sion of the Universe. In order to derive reliable absolute ages 
estimates, the most i mport ant parameter to determine is the dis- 
tance (see iRenziniL Il99ll) . Errors on the distance propagates 
through the age determination procedure in such a way that a 
0.08 error in the distance modulus (DM) propagates as a 12% 
error in the age. Different ages (younger or older) mean dif- 
ferent epoch of formation of the Galaxy. Therefore, it can be 
easily seen how the uncertainty in GGCs ages estimates sig- 
nificantly affects models of galaxy formation and cosmological 
evolution. 

The recent availabilit y of the HIPPARCOS catalog of 
trigonometric parallaxes JPerrvman et all 1 19971) has greatly 
improved the knowledge of GGCs distances. A large num ber of 
papers on the subject have been published since then (e. g.lRei 



1997 



2002; 



Gratton et allll997tlPont et allll998tlvandenberg etal. 
Gratton et all 120031) . most of which present larger dis 



tances, and hence younger ages, with respect to pre vious works 
(with the exception, among the cited works, of IPont et ail 
1998, who recover the pre-Hipparcos results of 1992). 



Vincenzo 



Testae-mail: 



Send offprint requests to 
testa@mporzio . astro . it 

* Based on data collected at ESO-La Silla, Chile, (GTO 63.L-07 17) 
and from HST archival data (GO 8310) 



The method used consists of a fit of the cluster main 
sequence to the subdwarfs main-sequence of corresponding 
metallicity, once the reddening has been subtracted. This pro- 
cedure is, of course, prone to uncertainties, both random and 
systematic: on the subdwarf side errors on parallax measure- 
ments, Lutz-Kelker correction, metallicity; on the cluster side 
reddening, metallicity, narrowness and accuracy of the mean 
ridge line. A further effect is represented by the presence of 
binary stars in the cluster CMD and in the subdwarfs sample. 
In the first case the determination of the mean ridge line is al- 
tered, if a proper fitting procedure is not taken into account. In 
the second, the magnitude of the field star is different from stars 
of analogous metallici ty. Th is issue has been w i dely discussed 
in, e.g ., lGratton et alJ Jl997l) : IPont et alJ Jl998l) : ICarretta et alJ 
(2000J), and will be analyzed in detail later. If HIPPARCOS 
gave the chance to significantly reduce the random errors on 
parallax measurements, nonetheless other sources of errors still 
remain. Metallicity, for example, represents a common prob- 
lem that significantly affects the determination of the distance 
modulus, because of either the spread in metallicity of the se- 
quence of local subdwarfs or the error on the absolute values. 
There are two strategies to go around this problem: to obtain a 
sufficiently large set of subdwarfs in the right metallicity range, 
or to apply a relation that transforms the colors of the stars of 
different metallicity. The second introduces a model-dependent 
relation and should be used with caution. Of course, the first so- 
lution is better, especially if the sample of subdwarfs draws a 
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main sequence down to very low masses. In fact, recent HST 
color-magnitude diagrams of GGCs reach very faint magni- 
tudes, hence small masses, but a similar sequence for subd- 
warfs is currently not available. iReid and Gizisldl998l) fit a dis- 
tance modulus to NGC 6397 by using a sequence of lower MS 
subdwarfs, but their work stands on the fact that the sequence 
of very metal-poor subdwarfs define a reasonably narrow se- 
quence. This is not true, instead, for subdwarfs of higher metal - 
licities (see their Fig. 2). On the other hand, the most widely 
used set of subdwarfs cover magnitudes ranging from My ~ 7 
and brighter, leaving uncovered the range between this value 
and the extreme subdwarfs. This is explained with the relative 
lack of halo stars near the Sun dReid and GizislH998l) . One of 
the possibilities we discuss in order to reduce the uncertain- 
ties is to perform a uniform study of a subdwarfs sample and 
a GGC by using the same telescope, instrumentation and, pos- 
sibly, observing nights, thus eliminating the systematics due 
to the relative calibrations among different instrumental se- 
tups. For this reason, we observed the target cluster and a sam- 
ple of subdwarfs in the proper metallicity range matching the 
cluster's with the same telescope+instrument (NTT+SUSI2) 
and then used archival data from HST to go faint and cali- 
brated the HST magnitudes into the ground system, obtain- 
ing a deep and fully homogenous data set. The selected ob- 
ject is the well known cluster NGC 5904 (M 5), for which a 
wide set of photometric an d spectroscopic studies is available 
in the l iterature. Mor eover, Sandqu^tetaljll99d S96), Reid 
Jl997h . lGratton et alJ dl997h IChabover et alJ (11993) and, more 
recentlvT Ivandenberg et alJ (120021) determined the DM of the 
cluster by using the trigonometric parallaxes of subdwarfs and 
various sets of models. These works will be often referred to in 
the following. 

NGC 5904 is a typical old halo clu ster of int ermediate 

metall icity -[Fe/H] 1.1 in the scale of Carretta and Grattoi] 

ill 9971) (CG), -1. 4 in the IZinn and WesJ Jl984 scale fZWI. 
A recent study of iRamirez and Cohenl J2003I) gives [Fe/H] 
—1.3. The morphology of the horizontal branch (HB) does not 
show signs of the second parameter effect and several age deter- 
minations are available. The structure of the paper is outlined 
as follows: Sect. 2 describes the observations and the archival 
material, Sect. 3 presents the reduction and calibration of the 
data, in Sect. 4 we present and analyze the subdwarfs sample. 
In Sect. 5 the cluster color-magnitude diagram (CMD) is pre- 
sented and discussed, and in Sect. 6 the DM-fit is analyzed. In 
Sect. 7 a final discussion is summarized. 

2. Data Set 

2.1. Observations 

Observations have been performed with the ESO-NTT tele- 
scope equipped with SUSI2 as a backup program during three 
different campaigns, in december 1999 (63.L-0717, subdwarfs 
+ cluster), december 2000 (66.D-0242, subdwarfs), february 
2002 (68.B-0061, cluster). The adopted strategy was to ob- 
serve both the subdwarfs sample and the target cluster with the 
same setup telescope+instrument+filters+detector and, possi- 
bly, within the same campaign, in order to remove any sys- 



tematics due to the use of different setups. The december 1999 
and february 2002 campaigns produced very good quality data, 
while the december 2000 set (two subdwarfs only) presented 
moderately poor photometric conditions. Seeing was not a ma- 
jor problem for observing the subdwarfs,because they had to 
be defocussed in order not to saturate the chip. This, coupled 
to the fact that the objects were observed in the central region 
of the frames, allowed to use exposure times long enough to 
have a negligible shutter effect. The typical observing time, per 
frame, was 0.5 s. 

The data set for the target cluster (NGC 5904 = M 5) was 
secured in december 1999, when the seeing was poor, together 
with the subdwarfs, and in february 2002 when seeing was con- 
siderably better, in order to have a comparison data-set. In both 
cases, the run was found to be photometric and, in fact, the run- 
to-run scatter in the magnitudes is almost absent, without any 
trend in the magnitudes and colors (see Fig|5}- The field was se- 
lected 2 arcmin North of the cluster center, to ensure an almost 
complete overlapping with the field observed with HST that 
have been retrieved from the archive (see below). Only images 
in V and I bands for the cluster were taken, because the space- 
based data-set consists of images in these two bands only. For 
the subdwarfs, also B frames were secured. 

2.2. Archival data 

A set of deep exposures was available from the STScI HST 
WFPC2 archive (GO 8310), in the two wide bands F555W (V) 
and F814W (I). A total of 6 F555W and 9 F814W frames were 
retrieved from the archive, and pre-reduced by using the HST 
pipeline, av ailable in STSDAS running under IRAF dTodvl 
19861 1 1 9931) and using the most recent calibration frames and 
tables (see the WFPC2 web page for details). 

3. Reductions and Calibrations 

3.1. Reductions 

Reductions have been done for all the data-sets by using IRAF 
and the digiphot packages daophot and photcal for building up 
the catalog and apply photometric calibrations respectively. In 
particular, since the telescope had been defocused when ob- 
serving the subdwarfs sample, for these stars aperture photom- 
etry has been used. For cluster images, both NTT and HST, we 
applied the standard daophot PSF fitting routines. In particular, 
the NTT data-set is characterized by a relatively high oversam- 
pling (FWHM ~ 5 pix), and in this case a gaussian model for 
the PSF plus look-up table of correction varying quadratically 
along the frames has been used. HST-WFPC2 frames, on the 
other hand, are undersampled (estimated FWHM ~ 1.6 pix), 
and in this case a moffatl5 function, with quadratic look-up ta- 
ble, gave better results. As a general strategy, object searching 
has been done independently on each frame, after registering 
all of them to a common reference frame (the best seeing one). 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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In the case of HST, all the images were almost perfectly regis- 
tered onto each other, and no shift or rotation was needed. We 
made use of the world coordinate system available in the FITS 
header in order to obtain later an astrometric calibration of the 
catalogs. 

3.2. Photometric calibrations 

For both data-set, aperture corrections were applied in or- 
der to report the PSF magnitudes to aperture magnitudes be- 
fore calibrating. HST-WFPC2 data were reported to the stan- 
dard aperture of .1 , f ollowing the calibration procedure of 
tHoltzmann et all 1 1995b . For ESO-NTT data, growth-curves 
(from DAOGROW) have been used and aperture corrections 
at infinity have been computed. For ground-based data, a set of 
Landolt standards has been secured for each run and calibration 
equations derived, in the form^ TableJQreports the calibration 
coefficients of all the runs. The equation for the V filter has 
been derived against B - V in the two first runs, as B,V,I data 
were available both for subdwarfs and for NGC 5904. In the 
last run, only NGC 5904 has been observed and only in V and 
I filters, hence we calibrated the V magnitude against (V-I). 

v = V + flo + a\ x (b - v) + a 2 x X v 

v = V + a + a[x (v - + a' 2 x Xy (1) 

b = B + bo + b\ x (b - v) + b 2 x X B 

i = I + cq + c\ X (v - + C2 X Xj 



Table 1. Calibration coefficients for the three runs. 



Filter 


Zero Point 


Color Term 


Extinction 




Dec. 1999 






B 


25.42(0.01) 


-0.05(0.01) 


0.23(0.01) 


V 


25.82(0.01) 


0.02(0.01) 


0.11(0.01) 


I 


24.74(0.01) 


-0.04(0.01) 


-0.04(0.01) 




Dec. 2000 






B 


25.86(0.03) 


-0.18(0.04) 


0.29(0.02) 


V 


25.94(0.04) 


0.00(0.04) 


0.19(0.03) 


I 


24.80(0.01) 


-0.17(0.01) 


-0.09(0.02) 




Feb. 2002 






V 


25.69(0.01) 


-0.01(0.01) 


0.15(0.01) 


I 


24.59(0.02) 


-0.03(0.01) 


-0.06(0.01) 



Figs.[2|2]and[3]show the difference between standard mags, 
and calibrated mags, of the Landolt stars used for calibrations 
in all the runs. As can be seen, the residuals (Mag Jfrf - Magy,-,) 
are on the average very small and with no residual trends with 
the colors. 

The HST-WFPC2 data- set has been cal i brated following 
the procedure describred in lHoltzmann et alJ j 1995b . After ap- 
plying the aperture corrections and c orrecting for g eometric 
distorsions and CTE effect, following Dolphin (2000), we ap- 
plied the published coefficients and calibrated the data-set both 
in the VEGAMAG and STMAG systems. While the V magni- 
tudes are almost identical in the two systems, the I magnitudes 



Figure 1. Residual between standard magnitudes and cali- 
brated magnitudes for the Landolt stars in the Dec. 1999 run. 
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Figure 2. Residual between standard magnitudes and cali- 
brated magnitudes for the Landolt stars in the Dec. 2000 run. 
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differ by a certain amount (around 1.2 mag, depending on the 
WF chip). 

The NTT data-set has then been used to report magnitudes 
in both system to the ground NTT system, in order to perform 
also an internal check. 

4. The subdwarfs sample 

The sample of subdwarfs used for this study has been se- 
lected from the HIP PAR COS l ist tak en from the papers of 
Gratton et and iReidl ill 9971) . The requirement was 

that the sample be as close as possible in metallicity to the tar- 
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Figure 3. Residual between standard magnitudes and cali- 
brated magnitudes for the Landolt stars in the Feb. 2002 run. 




0.5 

V-I 



get cluster, in order to reduce systematics introduced by any 
trasformation that would have been needed to report the values 
to the metallicity of M 5, and be distributed over a wide range 
of magnitudes and colors in order to overlap a sufficiently large 
portion of MS and ensure, in this way, a more robust fit to the 
cluster MS mean ridge line. The sample consists of 6 stars, 4 of 
which observed in the dec. 1999 run, the others in dec. 2000. 
Each star has at least three measurements, in order to reduce 
the error on the magnitudes and remove the effect of eventual 
bad pixels and/or fiat field dishomogeneities. 

Since these stars needed to be observed with very short ex- 
posure times, the shutter delay of SUSI2 has been taken into 
account. In particular, the typical exposure time was 0.5s, and 
the shutter delay for SUSI2, taken from the SUSI2 web page at 
ESO 2 . 

Table|2]shows the calibrated magnitudes and colors for the 
selected sample. The error indicated is the internal photomet- 
ric error, and is not representative of the total error. To obtain 
this last, uncertainties on the reddening and trigonometri c par- 
allax must be applied. As pointed out also by Gr atton et alJ 
J2003h . these stars have virtually zero reddening. The error on 
the trigonometric parallax, given in the literature, has been ap- 
plied following the propagation of errors. The Lutz-Kelker cor- 
rection has been co mputed on the basis of the discussion by 
ICarretta et alJ ll2000l) . by using the formula adopted by the au- 
thors (Eq. 1 in the cited work), because the sample had the same 
selection effect as the one reported in that paper. This value is of 
the order of 0.02 mag. The associated Malmquist bias, for the 
same reason, is negligible. By adding all the error sources, the 



average uncertainty increases to ~ 0.03mag but with a spread 
given by star-to-star differences. 

Since the selection has been done with the criterium on 
the metallicity explained above and requiring that the objects 
be, of course, visible from the site and in the period of obser- 
vation, only six subdwarfs have been selected, one of which 
saturated the I filter even with the shortest possible exposure. 
The effect of metallicity on the subdwarf fitting to the dis- 



cent papers (see,e.g 


, Gratton et al., 


1997; 


^ n 

Carretta et al., 


2000; 


Gratton et aill2003HReidlll997HPont et alJ 


1998) and will not 



be readdressed here. However, it is worth recalling that by se- 
lecting objects with a metallicity close to the cluster's one, the 
uncertainities due to the metallicity spread are reduced, and the 
correction to a common metallicity values is small. The metal- 
licity of the cluster has been recently recalibrated in the CG97 
scale in [Fe/H] = -1.10, while ZW84, and subsequent updates, 
give an estimate of [Fe/H] = -1.40. The metallicit y scale of 
the su bdwarfs, instead, has been recalibrated bv lGratton et alJ 
(120031) and foun d to be, on average, . 1 3 +0.04 dex smaller than 
the one given in lCarretta et alJ J2000h . This offset is to be taken 
into account in the DM determination, and will be discussed 
below. 



Table 3. Extra subdwarfs from literature. 


Went. 


M v 


(V-I)" 


[Fe/H] 


Source 


HIP57939 


6.61±0.02 


0.88 


-1.33+0.07 


G03,H 


HIP57450 


5.59±0.25 


0.63 


-1.26+0.07 


G97,H 


HIP74235 


6.74±0.08 


0.81 


-1.38+0.07 


G03,H 


HIP79537 


6.84+0.02 


0.94 


-1.39+0.13 


G97,H 


HIP7459 


5.43+0.23 


0.679 


-1.2 


R01 


HIP73798 


5.90+0.26 


0.727 


-1.2 


R01 


HIP89215 


6.50+0.26 


0.86 


-1.2 


R01 


HIP100568 


5.44+0.12 


0.678 


-1.00+0.07 


R01,G03 



" G97:Gratton et al. 97; H: HIPPARCOS catalog, for V-I color; 
R01: Reid et al. 2001; G03: Gratton et al. 03. 



5. The Color-Magnitude Diagram 

5.1. The NTT CMD 

The CMD derived from ground-based data (feb. 02) is shown 
in Fig.0] The diagram extends from the upper part of the RGB 
down to about V ~ 24. A reasonably large magnitude range on 
the MS is crucial in order to report the HST magnitudes (both 
VEGAMAG and STMAG) to the ground NTT system. The CMD 
from the dec. 99 run has been obtained under relatively poor 
seeing conditions and hence came out to be not suited to per- 
2 http://www.ls.eso.org/Telescopes/NEWNTT/refdata/param.ccd.html form an accurate calibration of the HST data to the ground sys- 
shows that the exposure is uniform down to 0.3s. Since the stars were tem. For this reason we adopted the feb. 02 data-set after check- 



placed close to the center of the chip, the resulting shutter effect turns 
out to be negligible 



ing that the two sets of measurements were self-consistent. 
Fig |5] shows the results of the check. 
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Table 2. The subdwarfs sample 



Went. 


B 


o-(B) 


V 


o-(V) 


I 


o-(I) 


[Fe/H] 


M v 




Run 


HIP70681 


9.881 


0.001 


9.302 


0.001 


8.577 


0.001 


-1.09+0.07 


5.72±0.16 


19.16 


Dec99 


HIP74234 


10.272 


0.001 


9.422 


0.001 


8.446 


0.001 


-1.28±0.07 


7.08±0.11 


33.68 


Dec99 


HIP74235 


9.821 


0.001 


9.059 


0.001 






-1.30±0.07 


6.74±0.08 


34.14 


Dec99 


HIP81170 


10.344 


0.001 


9.611 


0.001 


8.723 


0.001 


-1.14±0.07 


6.18+0.15 


20.71 


Dec99 


HIP3026 


9.722 


0.001 


9.269 


0.001 


8.710 


0.001 


-1.50 


4.15 


9.57 


DecOO 


HIP24316 


9.924 


0.001 


9.448 


0.001 


8.838 


0.001 


-1.44+0.07 


5.28±0.15 


14.55 


DecOO 





5.2. The HSTCMD 

The HST measurements have been averaged, recording also in- 
formation on the number of frames in which each object was 
detected in each filter. As final catalog we consider all the stars 



appearing in at least 80% of the frames (at least 5 V and 7 I). 
This choice is a good compromise between taking all the mea- 
surements and taking only the stars detected in all the frames. 
The first choice would contaminate the CMD with all the spu- 
rious detection that appear by chance on only one frame. The 
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second leads to missing faint objects, like WDs. We remark the 
fact that the brighter part of the catalog (down to V ~ 26.5) 
is not altered by this selection criterion, as expected. The HST 
CMD is reported in Fig.|6j where we show only the VEGAMAG 
magnitudes. 

The whole M5 CMD is reported in Fig.0 The HST data- 
set extends approximately from V = 19.5 to V = 28, but the 
top of the sequence is spread over because the chip is close 
to the saturation and will not be considered in the following 
analysis. F555W and F814W magnitudes have been reported 
to the ground based CMD via a relative calibration equation of 
the form (for the V magnitude, but it applies to the I magnitudes 
as well) 

V ground = VhST + Zp + Ct X (V ground ~ I ground) 

The output for the V magnitude turned out to be identical for 
STMAG as well as for VEGAMAG, while for the I filter we 
derived, as expected, different transformations. The two HST 
magnitudes systems have been checked for self-consistency af- 
ter reporting them to the ground NTT system. The two trans- 
formed magnitudes still show some residuals that, over a range 
20 < V < 23 is ~0.01 mag, with a shape that resembles that of a 
CMD. This is an indication that the color effects were not com- 
pletely eliminated with the transformations, probably because 
of the larger scatter in the ground-based CMD and the rela- 
tively small magnitude overlap of the ground- and space-based 
data. However, even at the fainter magnitudes, the overall dif- 
ference between the two magnitudes is always below 0.05 mag 
(at V ~ 27.5). In the following, we will adopt as HST trans- 
formed magnitudes the ones coming from VEGAMAG, because 
this system is closer to the ground Vega system, so that the re- 
quired transformation is smaller and thus less prone to residual 
uncertainties 3 . 



[S]shows the ridge line determined as described above, the fidu- 
cial by S9 6 and the ones derived with the same method by using 
the data o fljohnston and Boltd Jl998l) and R00, that are in close 
agreement with ours. The irregularities in the 



Figure 8. Comparison of ridge lines for M5. 
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ridge lines obtained with literature data, that can be seen at 
the base of the RGB, are due to the small size of the magnitude 
bins chosen to draw the fiducial. Since we are not interested 
in a careful analysis of those ridge lines, but only show them 
for sake of comparison, we did not apply any smoothing or 
more refined technique in order to improve it. It is, however, 
worth noting that our fiducial line extends about 5 magnitudes 
below previous works. This fact is crucial because, at least in 
principle, allows us to check model predictions down to the 
bottom of the MS, as well as to use very faint subdwarfs to 
constrain the distance fitting, if any (see below). 

5.4. Comparison with previous photometry 

Beside the comparison of the ridge lines, a further check has 
been done by comparing the NT T CMD with re cent literature 
results, namely R00, Johnst on and Boltd Jl998l) and S96. Fig. 
1101 shows the positioning of fields observed by us, the HST 
pointing, and the three comparison works. The plots in Figs. II II 
and El snow the comparison between our magnitudes and the 
cited ones, where an overlapping of the fields permitted a direct 
match of the objects. Our data show substantially a good agree- 
ment with S96 and R00 but there are a few issues that should 
be discussed: i) the average scatter is quite large, meaning that 
a trend with the position could be present. Actually, there is a 
small residual trend with position, with amplitude of the order 
of 0.02 mag. In the comparison with the two literature pho- 
tometries cited above the trend goes toward opposite directions, 
3 see also the "HST Data Handbook for WFPC2" at in the sense that the trend along the X-coordinate, in the refer- 
http://www. stsci .edu/hst/HST_overview/documents/datahandbnici; system of S96, resembles the trend along the Y-coordinate 



5.3. Mean Ridge Line 

The mean ridge line has been drawn on the faint part by us- 
ing HST data, on the upper MS until the sub-giant branch with 
the NTT plus lit erature data, on the upp er part of the RGB and 
the HB with the ( Rosenber g et alll200ol R00) data, since their 
RGB is very narrow and well populated. In general, in order to 
build the ridge line, the CMD has been 'sliced' in magnitude 
bins and, for each bin, the median point of the distribution has 
been determined, together with the associated scatter. Only the 
objects with the lowest errors have been selected, in order to en- 
sure that the uncertainty on the ridge be small. The zones with 
a strong curvature, i.e. at the bottom of the SGB, were treated 
differently, by applying a polynomial fit to that range and pro- 
jecting the stars along the curved abscissa, then re-fitting iter- 
atively the polynom. In this way, the CMD is sliced along the 
curve and it is possible to find the median point with high ac- 
curacy also in the most bent magnitude bins. Fig.|5]shows the 
ridge line of the lower MS with its uncertainty strip. The nar- 
rowness of the uncertainty strip will be used later to constraint 
the fit to the subdwarfs sequence. As a comparison term, Fig. 
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Figure 9. The ridge line of the lower MS of M 5 with its uncer- Figure 11. Comparison with the data-set of S96 
tainty strip. 
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Figure 12. Comparison with the data-set of R00 



Figure 10. Map of the fields covered by this work and recent 
literature. 
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in the R00 data-set, always in the reference system of S96 (all 
the data have been reported to the system of S96), and vicev- 
ersa. The amplitude of the trend is larger with S96 than with 
R00, and is absent if we compare the NTT 1999 data with both 
S96 and R00. The only valid explanation for this is that a resid- 
ual flat-field is still present in some data-sets with respect to 
the others. In fact, in the reduction of the Feb. 2002 data-set, 

we applied a more refined technique for flat-fielding the data, 

I . i j i 

described in Hainaut ( 1998), which takes into account also the 
lowest spatial frenquencies. The 1999 data-set has not been cor- 
rected for this effect. On the average, this removed a residual 
trend in the fluxes. Since this procedure has been applied with 
success in other works (Monelli et al. 2004, in preparation), we 




find no reason to discard the procedure or re-calibrate our data 
in the system of S96 or R00. However, the effect is visible es- 
pecially at fainter magnitudes (V > 19), while is almost absent 
on the RGB. 



6. Distance via subdwarfs fitting 

6.1. Error sources 

Before proceeding with the DM fitting, it is appropriate to dis- 
cuss the contribution of t he various error sou rces to the to- 
tal uncertainty on the DM. lGratton et al. I d2003l) report in their 
Table 1 a detailed list of all the error budget contributors before 
and after the re-analysis done in that paper. Taking that as a 
guideline, we can now see how the various uncertainties com- 
bine in our case. In Sect. 4 we have already mentioned some 
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Figure 13. Color-metal licity relation for the V-I color, taken 
from the models of Stra niero et al.Nl997h . 




of the typical uncertainties in the subdwarfs analysis, and the 
fact that, since our sample h ave be e n drawn from the same lists 
publis hed bv lGratton et alJ Jl997l):|Carretta et alJ J2000h : lReidl 
ill 9971) and, partially. IGratton et all J2003L we expect that the 
uncertainty given by the Lutz-Kelker effect, the Malmquist bias 
and, clearly, the parallax, are esse ntially the same. Hen ce we 
adopt the values given in Table 1 of lGratton etaflfe003h . 

Photometric errors are very small both for the subdwarfs 
(see Table |2jl and the cluster ridge line. For this last, in par- 
ticular, the very precise HST photometry coupled with a good 
ground-based calibration allowed us to obtain an average error 
on the ridge line of the order of 0.005 mag at (V - 1) = 0.7. The 
error from the calibration relation is also very small and of the 
order of 0.02 mag, since the photometric conditions of the cali- 
bration night were close to excellent. In any case subdwarfs and 
cluster have been observed in the most homogeneous possible 
conditions and thus, the errors on the photometric calibration 
are expected to be very small. A residual systematic trend could 
be present if the calibrated magnitudes were not reported to the 
standard system, but the obtained values are compliant with all 
the existing results, and thus we discard this possibility. 

Three other sources of error have to be now examined: the 
metallicity scale, the cluster's reddening, the presence of bi- 
naries. Their contribution affects both random and systematic 



6.1.1. Metallicity 



The metallicity scale is a major issue anc 
investigated in previous works (see e.g. G 


has been widely 
rattan et all I2003I 


Carretta et al., 2000; 


Reidlll997llPont et al. 


, 1998). There are 



two main aspects to take into account: i) the spread in metal- 
licity given by the distribution in [Fe/H] of field stars; ii) the 
relative metallicity scale used. In order to reduce the first ef- 
fect, the usual procedure is to apply a relation that transforms 



tha values to a common metallicity (the cluster's). In our case 
we decided to use only stars having metallicity values as close 
as possible to M 5, so that the transformation to a common 
metallicity is small. In order to estimate it, we used the re- 
lation betwee n [Fe/Hl and (V - I) at My = +6, taken from 
the models of Stra niero et al. Jl997l) . that is shown in Fig. 1131. 
A third order polynomial has been found to fit adequately the 
(V - I) vs [Fe/H]) relation in the range -2.3 < [Fe/H] < -0.5. 
The second effect stands on the systematic difference between 
the two mos t com monly used scales, the ZW and the CG. 
IGratton et alJ i 19971) use an homogeneous metallicity system, 
but o ur subdwarf s sample has been taken from that work and 
from Reid ( 1997), who used the ZW scale. Since our cluster 
has metallicity estimates [Fe/H] = -1.1 in the CG scale, and 
[Fe/ H] - — 1 .4 i n the ZW scale, we found it a ppropriate to re- 
port iReid ll 19971) values to IGratton et alJ 1 1997 ) i.e. fr om ZW 
to CG, and then selecting stars from the Reic ill 9971) sample 
with the reported metallicity close to the cluster's one. 

It should be noted that some o f the stars in the prese nt sam- 
ple have not been recalibrated by Gratton et al. ( 2003), hence 
a residual systematic spread might be present which affects the 
DM fitting as well, due to the co-existence of two slightly dif- 
ferent recalibrations. If the value of 0.13 dex is taken as aver- 
age recalibration offset, we estimate a residual systematic er- 
ror due to metallicity of +0.008 (at 1 cr) in the (V - I). Given 
the slope of the MS between M v = +5 and M v = +6, this 
propagates into an uncertainty of 0.04 mag. The uncertain- 
ly on the metallicity (usually 0. 1 dex) is, instead, added to the 
random error budget, and corresponds, following Gr atton et alJ 
(2003), to +0.08 mag. 



6.1.2. Reddening 

The cluster's reddening also plays a double role, since the asso- 
ciated error affects the global random error and the uncertainty 
on the absolute amount of reddening is a source of systematic 
error. We shall discuss quantitatively this point below when we 
present the fitting procedure. 



6.1 .3. Binaries 

The presence of binary stars affects both the subdwarfs and the 
cluster sample. Although we selected bona fide single stars, 
a residual quantity of binaries could still affect the subdwarf 
sample. We rely always on the f act that the sample ha s be en se- 
lected from the published lists of Gratt on et al.ldl997l) andlReid 
(I1997I) . and thus assume their figure for the contribution of sub- 
dwarf binaries to the total error, i.e. +0.02 (see IGratton et all 
2003). The effect of binaries on the cluster ridge line is slightly 
more complicated. The effect of the presence of binaries in the 
cluster is to spread the main sequence toward redder colors and, 
in extreme cases, to draw a separate seq uence, as has been de - 
tected in some open clusters (see, e.g.. iMarconi et all Il997l) . 
In our case, the HST main sequence is very narrow, and a well 
defined clear binary sequence would be visible. Since this is 
not the case, we tried to estimate the effect of a MS spread due 
to an unresolved sequence of binaries. In order to do this, we 
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"rectified" the main sequence and built a color histogram of 
the brightest portion of the HST MS, that has a uniform width 
all over the considered range (20 < V < 22). A gaussian pro- 
file has been fit to this profile, considering only the blue tail of 
the distribution, that is in principle unaffected by the presence 
of binaries. The fitted profile has then been subtracted from 
the original distribution. The residual red tail is an indication 
of the binary contribution. This procedure is quite robust and 
gives reliable figures for the uncertainties on the position of 
the MS locus in the CMD. From this, we derive an estimate of 
a possible systematic deviation of the MS ridge line from the 
"true" single star MS given by binaries of cr(V - 1) = 0.005, in 
addition to the random error due to the calculated spread of the 
MS. Again, this value in color propagates into a +0.025 mag of 
systematic uncertainty. Table0]below summarizes all the error 
sources involved in the DM fitting process. 

6.2. Distance Modulus Fit 

By taking into account all the considerations and the caveats 
described above, we used the subdwarfs sample to fit a DM 
to the cluster data. In order to do this, we adop t a value 
of red dening of E(B - V) = 0.035(±0.005) dCarretta et all 
2000, and refs. therein), corresp onding to E(V - I) - 0.056 , 
[obtained from the relation in iRieke and Leb ofskv (1985): 
E(V-I)= 1.6xE(B-V)]. The quoted uncertainty has 
been taken from the literature, and propagates into 
a 0.025 uncertainty in magnitude. The DIRBE maps 
techlegel et allll99 8) report a reddening value in the direc- 
tion of M 5 of E(B - V) = 0.036. cons istent with the absolute 
value given in ICarretta et alJ J2000). As shown in Fig. G3 
the brightest calibrating subdwarfs (HIP3026 and HIP24316) 
turn out to be too blue to fit directly the cluster ridge line, be- 
cause the cluster MS has already evolved toward the turn-off 
point. However, they are useful to draw the subdwarfs main se- 
quence locus to be compared with our data. Of the other three 
stars, HIP8 1 170 is known to be a spectroscopic binary, and this 
might justify its redder color. The other two, HIP70681 and 
HIP74234, have the smallest error and have then been used 
to fit the cluster ridge line. We fit the cluster fiducial line to 
the subdwarf sequence by shifting the fiducial according to the 
reddening (in (V— I)!) and extiction -3.1 xE(B - V)-, then with 
a least square fitting of the difference between the used subd- 
warfs and the cluster fiducial, with their errors. Fig.ll4lshows 
the results of the best fit, that turns out to be //o = 14.44. If 
a value of E(B — V) = 0.03 is adopted, which corresponds to 
E(V - 1) = 0.048, yu = 14.41 is obtained. The error on the 
fit obtained from the combination of the random error sources 
only is, in both cases, +0.09, the uncertainty in the metal- 
licity being the strongest contributor. In order to increase 
the number of subdwarfs suitable to form a fiducial reference 
sequence, we added to our sample other subdwarfs listed in 
iGratton et al] ( 1 1997b and lReid et aflfeOOll) . having [Fe/H] sim- 
ilar to that of M 5 and discarding known double stars. The 
added stars are listed in Tableland shown with different sym- 
bols (filled circles) in Fig. [21 It can be seen that most of the 
added stars draw a good reference main sequence, with the ex- 



ception of HIP57450 and HIP74235, which are considerably 
bluer. On the (V, B - V) CMD, however, they lie on the MS, 
and thus we are inclined to interpret their bluer colors as un- 
certainties in the I magnitudes, that has been taken from the 
original HIPPARCOS catalog. For this reason, the two outly- 
ing stars have been excluded from the fit. 

It can be seen that the added stars define three well sepa- 
rated sequences corresponding to the three metallicity ranges 
selected, from the lowest (triangles, bluest), to the intermediate 
(circles), to the highest (squares, reddest). 

With the extended sample of subdwarfs, our best estimate 
for the (dereddened) DM is no = 14.44 if E(B - V) = 0.035 is 
ado pted. If, in stead, E(B - V) = 0.03 is adopted, as reported in 
Idratton et ail lll997t) . the DM turns out to be /j. = 14.41, thus 
confirming the value obtained with the original, and homoge- 
neous, sample. 

The error associated to the fit above has been estimated 
from all the error sources discussed above and reported in 
Table |4] and considering random and systematic errors sepa- 
rately. From a statistical error propagation all over the involved 
quantites, the uncertainty associated to the fit turns out to be 
±0.09 mag. A systematic uncertainty, combined from redden- 
ing and metallicity scale, of ±0.07 mag must then be added. 
The global uncertainty is slightly larger than the one re- 
ported by Gratton et al. (2003), probably because the V - 1 
color is more sensitive to metallicity variation and thus the ef- 
fect of variation in [Fe/H] and E(B-V) are amplified with re- 
spect to the B-V color. 

7. Discussion and Conclusions 

This study was aimed at estimating the DM of M 5 through the 
main-sequence fitting method, using as calibrating sequence 
a set of local subdwarfs for which the trigonometric parallax 
from HIPPARCOS was known. In order to achieve our goal 
we used the same setup telescope+instrument+filters both for 
the cluster and subdwarfs data, namely the ESO-NTT equipped 
with the camera SUSI2 and standards ESO V and I filters. 

Later, a set of subdwarfs from literature has been added 
to the original sample in order to improve the sample size, 
and hence increase the accuracy of the fitting. All the com- 
parisons with recent data from literature show that our data 
are fully consistent with previous determinations, whose ac- 
curacy is known, both for the cluster and the local subdwarfs. 
A final sample of ~10 subdwarfs has been used to fit a DM to 
the cluster. By taking into account the two most commonly ac- 
cepted values for the reddening (i.e. E(B-V)=0.03 and 0.035), 
we obtained two very close values for the (dereddened) DM: 
HO = 14.41 ± 0.09 ± 0.07 and no = 14.44 ± 0.09 ± 0.07, re- 
spectively, where the first number comes from random error 
sources, the second from systematic e rror source s. These val- 
ues are in good agreement both with Reid ( 1997) (/^o = 14.45 
obtained with the subd warf s havi ng the smallest error on the 
parallax) and ICarrettae t al. (2000) (jx v = 14.59 correspond- 
ing t o uo = 14.48), with a slight preference for the shorter 
value. Vand enberg etaD 12002). instead, found fio = 14.38 with 
the use of non-diffusive isochrones and /jo = 14.30 with dif- 
fusive isochrones. In their work, they use a reddening value 
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Table 4. Contribution to the total uncertainty, in magnitudes, of the various error sources. Values are at lcr level. 





Field Stars 


Cluster 


Photometry 


+0.001 


±0.005 


Calibration 


±0.02 


±0.02 


Lutz-Kelker 


+0.02 




Malmquist 


Negligible 





Parallax 


±0.01 




Metallicity (rand.) 


±0.08 


±0.08 


Metallicity (syst.) 


±0.04 


±0.04 


Reddening (rand.) 


Negligible 


±0.025 


Reddening (syst.) 


Negligible 


±0.03 


Binaries 


±0.02 


±0.025 



Figure 14. Fit of the cluster fiducial to the subdwarfs sequence. 
The stars with the label on the right belong to the observed 
sample. The others have been taken from published works. 
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E(B - V) = 0.038 and a metallicity [Fe/H] = -1.4 for the clus- 
ter, discussing their preference for the ZW scale over the CG 
on the argument that the metallicity for M 5 is inconsistent with 
the metal determination for t he field stars, if the value is as high 
as -1.1. Moreover, Vandenb erg et alJ ( 120021) choose to not re- 
port the field stars to a common metallicity, but this effect acts 
on the spread of the relation rather than on the absolute value. 
In order to try to obtain a further constraint on the DM, 

we tried to use a set of extreme subdwarfs, using an analo- 

■ ■ ■ ■ 

gous approach to Reid and Gizis ( 1998), that has been found 
to be unapplicable because both the main sequence locus at 
the faint magnitude level of the extreme subdwarfs (My ~ 10) 
is dispersed and the subdwarfs of intermediate metallicity de- 
fine a sequence dispersed as well. In order to reduce the uncer- 



tainty, it would have been desirable to have a set of subdwarfs 
of fainter magnitudes (i.e. below My ~ 7) to take full advan- 
tage of the narrowness of the HST main sequence, but they lack 
in the HIPPARCOS sample. Moreover, at such faint magnitude 
level, the lower main sequence stars are very cool and the de- 
termination of their metallicity much more uncertain because it 
is very difficult to take into account all the molecular species in 
the athmosphere at these low temperatures. 
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